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Delay Minimizing User Association in Cellular 
Networks via Hierarchically Well-Separated Trees 
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Abstract —We study downlink delay minimization within the 
context of cellular user association policies that map mobile 
users to base stations. We note the delay minimum user asso¬ 
ciation problem fits within a broader class of network utility 
maximization and can be posed as a non-convex quadratic 
program. This non-convexity motivates a split quadratic objective 
function that captures the original problem’s inherent tradeoff: 
association with a station that provides the highest signal-to- 
interference-plus-noise ratio (SINR) vs. a station that is least 
congested. We find the split-term formulation is amenable to 
linearization by embedding the base stations in a hierarchically 
well-separated tree (HST), which offers a linear approximation 
with constant distortion. We provide a numerical comparison 
of several problem formulations and find that with appropriate 
optimization parameter selection, the quadratic reformulation 
produces association policies with sum delays that are close 
to that of the original network utility maximization. We also 
comment on the more difficult problem when idle base stations 
(those without associated users) are deactivated. 

Index Terms —Cellular network, user association, delay min¬ 
imization, quadratic program, linear approximation, hierarchi¬ 
cally well-separated trees. 

I. Introduction 


A. Motivation 

Kleinberg and Tardos JT'| investigated linearization of 
quadratic terms in the context of the metric labeling problem. 
Given a set of objects P and a set of labels L with pairwise 
relationships defined among the elements of both sets, metric 
labeling assigns a label to each object by minimizing a cost 
function involving both separation and assignment costs. Sepa¬ 
ration costs penalize assigning loosely related labels to closely 
related objects, while assignment costs penalize labeling an 
object with an unrelated label. In this paper, we seek to relate 
the metric labeling problem to the user association problem in 
a wireless cellular network, where association costs between 
an object (mobile user, MU) and label (base station, BS) are 
inversely proportional to achievable data rates between the 
object-label pair, while separation costs between labels can 
be thought of as penalizing traffic flows across backhaul links 
connecting BSs. In the approximation algorithm of Kleinberg 
and Tardos, the label distance metric is first embedded into 
a hierarchically well-separated tree (HST), which simplifies 
separation cost estimation. Embedding is commonly used for 
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tackling intractable combinatorial problems involving geo¬ 
metrical data |3j. Approximating the solution in such cases 
can be done in polynomial time once data is embedded into 
tree metrics. However, such embeddings tend to introduce 
distortion in most cases. Bartal @ introduced the notion of 
HSTs and proved the lower bound of distortion of embedding 
arbitrary metrics into HSTs’ to be O(logn) where n is the 
number of nodes in the source graph. Fakcharoenphol et al. 
|4j later introduced a deterministic algorithm for embedding 
arbitrary graphs into HSTs’ with a tight bound on distortion. 


B. Related Work 

The transition from traditional cellular networks to hetero¬ 
geneous networks (HetNets) has opened up many research and 
design questions including user association, these are gathered 
and detailed by Andrews |5] and Ghosh et al. |6]. Of interest 
in this paper are user association problem formulations, their 
complexity, and approximability. 

The objective of user association is typically maximization 
of user rates. The core problem often involves combinatorial 
optimization by mapping MUs to BSs 0-0- Methods to 
distribute, approximate, or heuristically solve are often the 
key differentiating factor, while several key model assump¬ 
tions also serve to differentiate approaches to the problem. 
Other approaches to cell association include using stochastic 
geometry to characterize the outage probability of distribution 
of rates of a typical user in the network IIM3- 

Fairness of user association schemes has been addressed by 
several papers 0, 0-0. Sang et al. 0 propose a cross¬ 
layer, scheduling and load-aware algorithm to maximize the 
network’s sum, weighted, a-proportional fair utility. Son et al. 


1131 propose off- and on-line algorithms to compute handoff 
and association rules to achieve network-wide proportial rate 
fairness across MUs with log-based utility. Bejerano et al. 
0 and Sun et al. p~6) explore association rules that promote 
max-min rate fairness across MUs. Kim et al. 0 propose 
a a-optimal user association rule related to a-proportional 
fairness ED that provides a tradeoff between individual user 
rate maximization and load balancing across BSs. Our work 
focuses on delay minimization, which is equivalent to utility 
maximization under a = 2 proportional fairness. 

User association policies may be either centralized or decen¬ 
tralized. Many centralized problem formulations are NP-hard 
and require simplifying assumptions or techniques. Corroy 


et al. 1181 pose a relaxation of a centralized user associa¬ 
tion problem into one that is quasi-concave and provide an 
upper bound on the optimal value of the original sum-rate 




maximization. Kim et al. eg propose a distributed iterative 
user association scheme whereby individual MUs and BSs 
take turns making association decisions (MUs) and advertising 
loads (BSs) and prove convergence to the optimal solution 
of the corresponding centralized optimization problem. Shen 
et al. |9| develop a distributed pricing-based association al¬ 
gorithm that is based on the technique of using coordinate 
descent method on the dual of the original utility maximization 
problem. Ye, et al. © propose a distributed user association 
algorithm based on primal-dual decomposition of an initially 
centralized network utility maximization problem. 


C. Contributions 

In sjn] we introduce our downlink cellular network model 


and pose a sum-rate maximization problem (|8j. In f Ill 


we 

pose an alternate sum-delay minimization problem m and 
show that it is quadratic but non-convex. In §IV| we propose 
a split-term delay minimization problem ( fl3| that captures 
user rate maximization with BS congestion minimization. We 
show that the quadratic congestion term is amenable to linear 
approximation via HST embedding (Prop. [3} and bound the 
distortion that this technique introduces into our problem 
formulation (Lem. [TJ. Thus, a linear approximation ( [T6| is 
proposed in place of the split-term quadratic delay mini¬ 
mization formulation. In fjv] we discuss permitting idle BSs, 
those without any associated MUs, to deactivate and reduce 
interference. In f VI we provide a numerical comparison of 
the original combinatorial user association problem, related 
quadratic relaxation, split-term reformulation, and linear ap¬ 
proximation. Finally, VI1| concludes our work. 


II. Downlink Rate Maximization 

Consider a set of BSs B at locations {y a \a £ B} and a 
set of MUs U at locations {y p \p £ U}, where both sets exist 
within a bounded arena A C M 2 . A user association policy 
is a mapping / : U —► B that assigns each MU to exactly 
one BS. Although modem cellular standards allow MUs to 
associate with multiple BSs simultaneously, e.g., coordinated 
multipoint (CoMP), in this paper we ignore this generalization. 

We model the instantaneous downlink rate of MU p from BS 
a by the Shannon rate of their point-to-point channel, treating 
interference from other BSs as additive white Gaussian noise: 


r p ,a = log(l + sinr po ), 


(1) 


with signal to interference plus noise ratio (SINR): 


. = Pag(ya,y P ) 

pM Es >&B \aPb9(yb,y P ) + N ' 


( 2 ) 


determined by the BS transmission powers {p a \a £ B}. the 
channel attenuation function g(y,y') between locations y,y ', 
and the noise power N > 0. We model channel attenuation 
using large scale pathloss with exponent 7 > 2: 


g{y 1 y') = d{y 1 y') 1 = \\y-y'\\ 2 1 . (3) 


An assignment / \U —y B induces a partition on the set of 
MUs U , indexed by the set of BSs: ld a (f) = {p\f(p) = a,P £ 


U}. a £ B. Denote the resulting partition cell cardinalities 
as K a (f) = |if a (/)( and let n p (f) denote the cardinality 
of the cell to which MU p is associated under mapping /, 
i.e., K p (f) = n /(p)(/). Each BS is assumed to multiplex its 
resources (i.e., time for TDMA, frequency for FDMA) fairly 
across all of its associated users, so that the actual downlink 
rate to MU p from BS a under assignment / is: 


R P (f) = 


r P j( P ) 

«p(/) ’ 


( 4 ) 


This assumption of equal allotment of resources across MUs 
by each BS is only a convenience; future work will investigate 
expanding the model to incorporate potentially non-uniform 
allotments, as is done in |8J, although we note their model 
recovers uniform allotments as the optimal allotment (c.f. 
Prop. 1). The sum downlink rate of the network is: 


R(f) = Y, R p(f)- ( 5 ) 

p&A 


A reasonable objective of the downlink user association prob¬ 
lem might be to find the association / that maximizes the sum 
downlink rate of the network: 


maxf?(/). (6) 

Remark 1. The user association problem in (© intuitively 
trades off the desire for a MU to associate with its strongest 
BS, (a*(p ) £ argmax og g sinr p o ), vs. the least-loaded BS, 
(a*(p) £ argmin a6B K a {f)). 

We can express 0 directly as a combinatorial optimization 
problem using x = (x Pia ,p £ U, a £ B) as the vector of 
binary variables that together encode the assignment /. We 
set x ]Ka = 1 iff MU p is assigned to BS a. The set of feasible 
assignments under single-association is denoted X 7/ ;. 


: = < x : x £ {0,1} 


\U\\B\ 


i / J 'd'PiCL 
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( 7 ) 


We can write the rate of MU p as f p j( p )/n p (f) = 
Ea G eOW K, a (f))x p ^a- Next, we can represent the occupancy 
of BS a as rt a (f) = E q e_u x q,a- Substituting both expressions 
into (|6]» yields the following nonlinear integer optimization: 


max y 

xGAz 

p£U a(=:B 


Pi a 


E 


qeu x i, a 


h P:CL 7 


( 8 ) 


where the nonlinearity stems from BS resource multiplexing. 

As Ye et al. |8 | state, a brute force approach to the combina¬ 
torial user association problem has a complexity of 0(|6|l w l). 
Following their strategy, we work towards alternative problem 
formulations to ([8} that are easier to approximate or estimate 
for larger problem instances. 


III. Downlink Delay Minimization 

We now consider downlink delay minimization [ 191. which 
falls under a broader class of network utility maximization 
problems, particularly those that employ o-proportional fair¬ 
ness measures [|17)|, |2()|. When a = 2, each MU’s rate 







is assigned a utility equal to the negation of its reciprocal: 
— 1 /R p {f). We note that log-utility (a = 1) combined with 
fractional user association constraint relaxation, turns the net¬ 
work utility problem into one that is convex Q. We shall find 
that delay minimization results in a quadratic problem upon a 
similar relaxation, albeit a non-convex one. 

The sum downlink delay of the network and the sum 
downlink delay minimization problem are: 

D if) = E ~r~T 7\ m } n (9) 

p&A U P [J> 1 


We now provide an equivalent combinatorial representation 
of the delay minimization problem, analogous to ([8]) for the 
rate maximization problem: 


min 

xeAz 


EEE 

a£l3 p€U qEU 


-%v.a%a 


Vi a 


( 10 ) 


noting that it is quadratic in the assignment variables x. 

We relax the integrality constraints in (j9j) and let X- : denote 
the new feasible set of fractional assignments: 


T’ R =jx:xe [0,l]l w H B l,^a: p , o = l,Vpew|, (11) 

l aeB ) 


yielding a quadratic optimization problem over Ar: 

S EEE—W 9 ,a- 

e “ aeB pelt qeu rp,a 

We find that, in general, this problem is non-convex: 


( 12 ) 


Proposition 1 (Non-Convexity of One-Term Delay Minimiza¬ 
tion). Problem m is non-convex. 


Proof: Omitted for brevity. ■ 

To recap, we began with the natural but difficult to solve 
combinatorial optimization problem (equivalently, ©). 
then considered the modified a = 2 proportional fair combina¬ 
torial optimization problem (|9]i (equivalently, (| 1 ()[>), for which 
integer relaxation yields the non-convex nonlinear program 
( fl2| >. The similar agenda in using a = 1 and integer 
relaxation yielded a convex program, illustrating an important 
difference between a = 1 and a = 2. In <JTvl we consider a 


natural variant of (12 1 which we show to be a convex program. 


IV. Linearization of Min Delay via HST Embedding 


While we currently do not know how to directly convexify 
or linearize m we propose splitting the congestion term into 
separate terms; an assignment cost (first term) incurred by each 
MU and a congestion cost (second term) measured across pairs 
of MUs assigned to the same BS: 

min (1-/?) y'V'- x P , a +P V V Y' x p , a x q , a (13) 

lL — J r D a z ^ zz ' 

p£U a£t3 a£l3 p£U q£U 

where £ [0,1] controls the relative weighting between 
assignment and congestion costs. The cost of assigning MU p 
to a BS a is proportional to the instantaneous delay, 1 /r PiQ . 
Attempting to minimize assignment costs j3 = 0 would 


result in assigning each MU to the BS providing the lowest 
instantaneous delay. The cost of congestion is proportional to 
the number of MU-pairs associated with the same BS; i.e., 
for each MU pair (p, q), we pay a cost of J2aeB x P ,aX q , a - 
Assigning all MUs to one BS will make the sum 0(|W| 2 ), 
while attempting to minimize this congestion cost would result 
in an even distribution of MUs across BSs. Both cost terms 
maintain the tradeoff (Rem. [T]) between associating with the 
strongest BS vs. associating with the least congested BS. 

This formulation m is in part motivated by |8), where 
the choice of a logarithmic utility function naturally results 
in a similar split-term formulation. The split term formulation 
provides two benefits: i) the formulation becomes convex, and 
ii) the formulation is amenable to linear approximation via 
HST embedding. 

Proposition 2 (Convexity of Two-Term Delay Minimization). 
Problem m is convex. 

Proof: Omitted for brevity. ■ 

We now detail steps taken to restate m in linear form 
( |~i~6| ) by reformulating the quadratic congestion term. In its 
current form, each pair of MUs p and q assigned to BS a will 
contribute to the overall congestion sum. Instead, we define 
the congestion term using a unit distance complete graph G = 
(B. E ) whose vertices B correspond to BSs. Minimizing the 
congestion is then formulated as assigning MUs to BSs in 
G such that sum of the pairwise distances among MUs as 
measured in G is maximized. 


Proposition 3. A set of MUs IA can be uniformly distributed 
over a set of BSs B by solving: 

xfA’s 2 EEEE d'(a,b)x p , a x q , b (14) 

Z Z pGUqGUaGBbGB 

where d'(a, b) is the distance metric defined over graph G. 


Proof: Omitted for brevity. ■ 

Using an argument similar to 0. the quadratic term ( p~4| ) 
can be linearized by using the embedding of G into an HST 
T as follows: 


max 

xeXz 


s.t. 


) E E lTXeT 

e=(p,q) TCT 
p,q&U 

| x pT - x qT \ = x eTl e = (p, q), TCT. 


(15) 


Here x p t is the probability of MU p being assigned to a BS 
located in subtree T, i.e., x p t = Y^aeUT) x p,a, where L(T) is 
the set of BSs included in subtree T. We also note that d 1 (a, b) 
in the formulation (14 1 is replaced by Y2t ( t\ x p t — x q r I- 
We can now restate ( | 1 3| using the minimization form of 
© as follows: 


min 

xGAi 


S.t. 


- 0) E E 7-**“ - 0 E E hx e T (16) 

p£U aeB rp ’ a e=(p,q) TCT 
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Relaxing the integrality conditions to x £ Tt* in ( [T6| , we 
obtain a linear program. Solving this linear program yields a 
fractional solution to the original problem ( fl3| l. The fractional 
solution can be rounded in the same manner as fractional 
assignments in the metric labeling problem [TJ- We also note 
that the embedding of BS graph into T introduces constant 
distortion in the computation of congestion. 


Lemma 1. The solution of linear program (761 has 0(1) 
distortion. 


Proof: Omitted for brevity. ■ 

In summary, the split-term delay minimization formulation 
<0 is convex (Prop. |2j. We also note that the linear approx¬ 
imation of ( | 131 via HST embedding may be performed with 
a constant distortion guarantee (Lem. [TJ. 

V. Deactivation of Idle Base Stations 


In the problem formulations discussed thus far, the downlink 
interference seen by MU p when associated with BS a is 
independent of the assignment /. That is, although MU p’s 
association via / determines which BS carries its signal, the 
interference from the other BSs is assumed fixed. Although 
this independence is reasonable when all BSs have a non¬ 
empty associated set of MUs ( n a (f ) > 0,Va £ B), it is 
unreasonable for a BS to be assumed to transmit energy as in¬ 
terference when its association set is empty. The occurrence of 
empty association sets may be highly likely in heterogeneous 
networks containing a large number of small femtocells. If we 
define £>(/) = {a : n a (f) > 0} C B as the set of active BSs 
under /, then the SINR (|2]i will depend upon / as: 


sinr p (/) 


_ Pf(p)9(yf(p)iUp ) _ 

E**(/)\,(p) Pbdivb, Up) + N 


(17) 


where the interference in the denominator now only comes 
from active BSs B(f). 

The instantaneous rate ([TJ becomes f p (f) = log(l + 
sinr p (/)), the downlink rate ([dji to p becomes R p {f) = 
f p (/)/ K /(p)(/)’ lh e sum downlink rate of the network (|5j 
becomes R(f) = Y^ P eu 7? P (/), an d the sum downlink delay 
of the network § becomes D(f) = J2 P& u V^p(/)- Finally, 
the downlink user association problems ([6j and ([9j when idle 
BSs don’t transmit are: 


max R(f), min D(f). (18) 

Remark 2. When idle BSs don’t transmit, the MU decision¬ 
making tradeoff (Rem. [Tjl gains another dimension, i.e., MUs 
also seek to avoid the ‘interference activation cost’ by minimiz¬ 
ing the number of active BSs. To be precise, MUs on one hand 
wish to load-balance themselves across BSs (i.e., selecting 
associations f that induce equitable partitions U(f)), while 
on the other hand they wish to aggregate so as to minimize 
interference (i.e., selecting associations f with maximally 
imbalanced partitions lA{f)). 

We introduce binary BS variables z that take value one iff 
one or more MUs are associated with the corresponding BS: 


h 
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Fig. 1. Linear network with deactivation: assignment fa : {(p, a), ( q , a)} on 
the left, assignment fa : {(p, a), (g, 6)} on the right. Desired signals in blue, 
interference signals in red. 


z a (x) = 1 {J2 P eu x p, a > 0},Va £ £>■ Th e instantaneous rate 
from a BS a to MU p, now a function of the association 
variables, can be written as: 


t>,a(x) = log 1 + 


p a g(ya,y P ) 


J2b£B(f)\a Pbdiyb, y P )z b {x) + N 


09) 

Both problems in ( |T~8j can be written as (|8j and (10 1 with r pa 
replaced with r pM ix) as shown in (] 1 9[>. 


Under the addition of variables z a (x), problem (12 1 is 
no longer quadratic and provides stronger motivation for 
reformulating the problem in a more digestible form. One 
might include a term that penalizes the activation of BSs that 
are close to one another (interference avoidance): 


EE f(pa,Pb)d(a,b)z a (x)z b (x), (20) 

ae B bGB 


where d(a, b) is the distance between BSs a and b and 
f(p a , Pb) is a function of the transmit powers of both BSs that 
does not depend on z or x. Note that this term is quadratic 
in variables z, but nonlinear in x. We leave the maturation of 
this formulation for future work. 

In general, the deactivation of BSs improves network rates 
and delays due to a decrease in interference. In the following 
small network example, we show that there exist scenarios in 
which the optimal max sum rate and min sum delay association 
policies do indeed benefit from BS deactivation. 


A. Example: Linear Netw’ork with Deactivation 


Consider a simple network consisting of two BSs B = 
{a, b} and two MUs U = {p, q}. Let the BSs be located at 
y a = — d and y b = d with unit power p a = p b = 1. Let the 
MUs be located at y p = —Sd and y q = Sd with S £ (0,1). The 
parameter <5 controls how far the MUs lie away from the origin. 
Fig. □ provides the network layout for two of the possible 
four assignments in this example. Observe that assignment 
/*3 : {{p, b), (q, o)} is inferior to fi\ both assign one MU to 
each BS, but f 2 universally minimizes MU to associated BS 
distances, thus increasing rate and decreasing delay for all 
MUs. Additionally, assignment fa : {( p , 6), (q, 6)} is identical 
to /i due to symmetry of the example. Thus, it suffices to 


optimize problems (181 under reduced assignment set {fa, f 2 }. 


Proposition 4 (Rate Maximization on Linear Network with 


Deactivation). The rate maximization problem in ( 18 1 on the 
linear network in Fig. [7] can be expressed as a threshold 
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Fig. 2. Linear network deacti vatio n example for 7 = 3 and d = 100 meters. 
Plotted are (to p-le ft) l.h.s. of HD vs. S for constant N = —90 dBW; (top- 
right) l.h.s. of \22) vs. N for constant 5 = 0.5; and (bottom) (<5*, N*) pairs 
on the decision boundaries of HD and (22). 


comparison between assignments fi and fa: 


A + d-g)~7 i/2f 1+ (i+g)-7 

' Ndi ) ^ ' Nd-y ) 
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Proof: Omitted for brevity. 
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1. 


( 21 ) 


Proposition 5 (Delay Minimization on Linear Network with 
Deactivation). The delay minimization problem in ( |7<S| ) on the 
linear network in Fig. [7] can be expressed as a threshold 
comparison between assignments fi and / 2 -‘ 






> /i 1 


'i , (1- 


Proof: Omitted for brevity. ■ 

In Fig. [2] we demonstrate parameter regimes that favor 
either assignment fi or assignment fi- In the top-left plot, 
we see that the closer the MUs are to the origin (<5 small), 
assignment fi becomes the optimal assignment. In the top- 
right plot, we see that for lower background noise (N small), 
assignment fi becomes the optimal assignment. In these 
cases, SINR maximization/interference reduction takes priority 
over load balancing. In fact, we see a partitioning of the 
(<5, N) parameter space based on which assignment is optimal. 
Interference from additional active BSs tends to have less of an 
impact when i) background noise is dominant (large N), or ii) 
distances from each MU to the two BSs are highly dissimilar 
(large 5). Note the latter produces high SINRs despite the 
activation of both BSs. In this regime, we see that load¬ 
balancing assignment /2 is favored over assignment fi. 


VI. Numerical Results 

We first provide a comparison of a subset of the problem 
formulations presented in this paper on a small scale network 
amenable to combinatorial evaluation. We consider a 100 m 



Fig. 3. Reported optimal user assignments, after rounding to an integer 
solution, on a constructed topology (base station as red circles, users as blue 
crosses). Shown are the combinatorial one-term sum rate maximization (ClR) 
H) (top-left), combinatorial one-term sum delay minimization (C 1 D) {To} (top- 
middle), quadratic one-term sum delay minimization (QlD) HD (top-right), 
quadratic two-term delay minimization (q2d) HD (middle row), and linear 
two-term delay minimization (l2d) \\6\ (bottom row). Congestion costs are 
controlled using (3 = {0,0.5, 0.99} from left to right. 


by 100 m area consisting of 4 MUs and 4 BSs. The positions 
of MUs and BSs are chosen to illustrate the effect of the 
parameter 3 on the split-term delay minimization formulation. 
We consider unit BS power and a pathloss constant of 7 = 3. 

For this network topology (Fig. [7}, we solve and com¬ 
pare several problem formulations (detailed in the caption of 
Fig- 0 - Due to the non-convexity of QlD, the reported solution 
may only be locally optimal. For each of the relaxed problem 
formulations (QlD, q2d, and l2d), we round the reported 
fractional solution into an integer solution by associating each 
MU with the BS for whom its fractional association was 
largest. 

Fig. [D displays the assignment solutions reported by each 
problem formulation for the constructed network topology. 
In this case, ClR (top-left) avoids congestion at the BSs 
completely. Alternately, clD (top-middle) and QlD (top- 
right) avoid the high-delays associated with longer MU-BS 
distances. Finally, q2d and l2d (middle and bottom rows, 
resp.) highlight the tradeoff between instantaneous delay term 
minimization /3 = 0 and congestion term minimization /? = 
0.99. Almost-pure congestion minimization (/) = 0.99) causes 
each BS to be loaded evenly with one MU each, while pure 
instantaneous delay minimization (f3 = 0) results in assigning 
MUs to BSs purely based on SINR maximization; all MUs are 
assigned to the single, closest BS regardless of congestion. 

An appropriate choice of /3 strikes a balance between both 
objectives and results in an association similar to that of 
ClD (top-middle). Fig. [4] compares the sum rates and delays 
achieved by each of the formulations shown in Fig. [3] as a 
function of 3. First, note that ClR, ClD, and QlD formulations 
are independent of 3- Next, oddly, q2d with nearly-pure 
congestion minimization (/3 —> 1.0) comes close to ClR as a 


















































2.5 



25 

CQ. 

t 20 


]' . L7 M . 

..R_ 


& 

^ 15 



sum 



5 



ClR - ClD - QlD .. Q2D . L2D 


in this case, the SINR from a BS to a MU is additionally a 
function of the association map. Future work includes exam¬ 
ination of the proper tradeoff between instantaneous MU rate 
maximization and BS congestion minimization (how to select 
/?), the dependence of user association problem complexity on 
proportional fairness parameter a, and further investigation of 
idle BS deactivation. 


Fig. 4. The resulting sum rate (left) and sum delay (right) for the constructed 
topology shown in Fig. [3] Legend corresponds to those in Fig. [5] 



Fig. 5. The resulting sum delay per MU (left) as a function of the number of 
BSs. Formulation Q2D is sampled at 0 values from 0 to 1 in increments of 
0.05, the 0 value from this set that minimizes the sum delay (left) is reported 
in the figure on the (right). Legends correspond to those in Fig. [5] with the 
addition of distributed heuristic assignment policies mindist and maxSINR 
which minimize MU-BS distances or maximize MU-BS SINRs, respectively. 


means of maximizing rate. Lastly, we focus on the sum delay 
(right) and see that for an appropriate choice of /3, q2d and 
l2d can match the performance of QlD, and all three come 
close to or meet the minimum sum delay reported by ClD. 

Finally, Fig.[5]shows the delay performance of the quadratic 
problem formulations as the size of a randomly generated 
network grows from 1 to 20 BSs while the MU to BS ratio is 
fixed at 5:1. BS transmit power values are assigned uniformly 
at random from the discrete set {50,125,250}. We see that 
QlD and q2d achieve lower sum delay (left) than commonly 
studied distributed heuristics mindist and maxSINR. Addition¬ 
ally, for the sampled values of (3 for q2d, we found that the 
delay-minimizing /3 (right) tended vary between [0,0.5] as the 
number of MUs and BSs increased. 

VII. Conclusion 

In this paper, we posed delay minimization via user asso¬ 
ciation as a quadratic network utility maximization program. 
While this quadratic representation is in general non-convex, 
we proposed an alternate split quadratic objective function 
that attempts to capture the inherent tradeoff: association with 
a BS that provides the highest SINR vs. a BS that is least 
congested. We were able to extend the technique of metric 
embedding via HSTs to our split-term quadratic program, 
in which the quadratic costs of BS congestion were linearly 
approximated with constant distortion. We provided a numer¬ 
ical comparison of several problem formulations and found 
that with appropriate optimization parameter selection (/?), 
the split-term quadratic formulation produced sum delays that 
were close to that of the original network utility maximization 
problem. We also commented on the more difficult problem 
when idle BSs (those without associated MUs) are deactivated; 
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